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Experimental observation of the unidirectional motion of a topological soliton driven by a bi- 
harmonic ac force of zero mean is reported. The observation is made by measuring the current- 
voltage characteristics for a fluxon trapped in an annular Josephson junction that was placed into a 
microwave field. The measured dependence of the fluxon mean velocity (rectified voltage) at zero dc 
bias versus the phase shift between the first and second harmonic of the driving force is in qualitative 
agreement with theoretical expectations. 
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The role played by solitons in various physical systems is important and commonly known. Experimental observation 
of dynamical effects produced by solitons is in many cases difficult because real systems are very often far from the 
idealized mathematical models which give rise to soliton solutions. Solitons of the most common type are topological 
kinks, a well-known example being a magnetic flux quantum (fluxon) in long Joscphson junctions (LJJs) 1,2 . A fluxon 
in LJJ can be easily driven by the bias current applied to the junction. The motion of a fluxon gives rise to the dc 
voltage V across the junction, which is proportional to the fluxon's mean velocity. Varying the dc bias current /, one 
can produce a dependence V(I), which is the main dynamical characteristic of the LJJ. An experimentally obtained 
V(I) curve allows to identify many interesting features in the dynamics of fluxons trapped in a LJJ 2 . 

A microwave field applied to a LJJ gives rise to an ac drive acting on the fluxon. In a spatially homogeneous 
lossy system, an ac drive may only induce an oscillatory motion of a kink, which is hard to observe in LJJ due to 
the absence of dc voltage. It was first predicted 3 and then observed experimentally 4 that an ac drive can support 
motion of a kink with a nonzero average velocity (v) in a system with a periodic spatial modulation. Thus, if an 
inhomogeneous spatial potential is present, the fluxon motion can be locked to the external ac drive frequency. In 
this case the vertical steps ( "Shapiro steps" ) appear in the dependence of the average fluxon velocity (v) on the dc 
bias. The fluxon velocity in this frequency-locked state can be either positive or negative, depending on the initial 
conditions. Another interesting effect for a homogeneous system, in which the ac-driven fluxon drift may occur in 
either direction but is not frequency locked, was recently reported by Goldobin et al. 5 . 

Transport phenomena in nonlinear systems can also be induced by the so called ratchet effect. This phenomenon has 
been studied in a large variety of physical systems 6 and manifests itself by an unidirectional motion under the influence 
of forces with zero mean. The underlying mechanism is the breaking of those spatial and/or temporal symmetries 7 
of the system, which connect trajectories with specular velocities (as well as the possibility of phasedocking of the 
particle dynamics by the external field 8 ). So far, great attention has been devoted to the study of the soliton motion in 
ratchet potentials. Josephson fluxon ratchets have been proposed and designed using either continuous long junctions 
of special shapes 9,10 or discrete junction arrays 11 ' 12 , including the quantum case 13 . In these structures, an ac drive 
applied to a fluxon placed into an asymmetric (ratchetdike) spatial potential leads to the rectified voltage at zero bias. 
The mechanism underlying soliton ratchets in asymmetric potentials and driven by harmonic forces was investigated 
for the double-sine-Gordon system which can be mapped to an array of 3-junction SQUIDs 14 . The possibility of a 
ratchet-like effect induced by more complex ac drives but in the absence of any spatial ratchet potential was also 
demonstrated 15, 16 . 

The importance of the ratchet-like effect induced by a bi-harmonic force in the context of Josephson junctions resides 
in the fact that it gives a non-zero voltage state in absence of dc bias. The resulting current-voltage characteristics 
of a uniform LJJ is fully controlled by the phase difference between the two microwave field components driving the 
junction, which is a nice feature for practical applications. In a more general context, the experimental observation of 
soliton ratchets induced by asymmetric forces, first reported in this paper, opens new perspective for soliton transport 
in other physical systems. In contrast to soliton ratchets induced by asymmetric potentials, which are generally more 
difficult to implement and to control (a change of the potential requires a structural change of the system), soliton 
ratchets induced by asymmetric external fields can be implemented in any physical system described by soliton 
equations. 

The aim of this Letter is to present the first experimental observation of the rectified dc voltage induced by the 
unidirectional motion of a fluxon driven by bi-harmonic microwaves of zero average, applied to a spatially uniform 
long Joscphson junction. Our results qualitatively confirm the theoretical predictions 15 ' 16 . 

Experiments have been performed using Nb/Al-AlO^/Nb Josephson junctions 17 . Wc have measured two samples 
which have shown similar behavior; here we present data for one of them having the mean diameter 2R = 100 /im 
and the annulus width W = 4 /an. The dc bias current was applied to the ring-shaped junction via the leads having 
the so-called Lyngby geometry 18 (schematically shown in Fig. 1). In this configuration the bias electrode width equals 
the ring diameter, and the bias current is distributed very uniformly around the sample. This is due to the fact that 
the current distribution carried by a superconducting strip is larger on the edges and its profile (R 2 — y 2 )^ 1 ^ 2 (where 
y is the transverse coordinate 19 ) exactly matches the uniform biasing of the junction. 

The critical current I c of this junction in the state with no trapped fluxons was found to be about 0.9 of the 
maximum expected value, corresponding to a nearly uniform current flow over the junction area. The junction has 
the critical current density of about 1 kA/cm 2 , which corresponds to the Josephson length of Aj w 12 /im and plasma 
frequency of about f p = 120 GHz. This implies the normalized junction length 2-kR/Xj = I m 26 and width W/Xj < 1, 
i.e. the junction can be regarded as long and quasi-one-dimensional. The measurements were done at the temperature 
4.2 K, using a shielded low-noise measurement setup. The radio-frequency (rf) current was supplied by means of an 
open-ended coaxial cable antenna placed above the junction. The antenna was oriented parallel to the dc bias current 
leads, therefore the inductively coupled rf current was flowing in parallel to dc current. The applied ac power was 
supplied by the two microwave sources (HP8672A and HP83620A, combined via a directional coupler) having the 
frequencies f\ = f = uj/(2-k) and f 2 = 2/ with / in the range between 1 and 5 GHz. The sources were phase locked 
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to a common 10 MHz reference signal. The phase shift 9 between the sources was controlled electronically via the 
phase of the reference signal and monitored using a network analyzer. The ac power levels P\ and P2, mentioned 
below, pertain to the input at the top of the cryostat. 

Trapping of a fluxon in the junction was achieved by cooling the sample below the critical temperature T c w 9.2 K 
for the transition of Nb into the superconductive state, with a small dc bias current applied to the junction 20 . The 
residual fluxon depinning current idep was found to be rather small, less than 3% of the Josephson critical current 
I c , measured without the trapped fluxon. As a fluxon can only be trapped by junction's local inhomogeneities in the 
absence of the magnetic field, this indicates a fairly high uniformity of the junction. 

In our experimental configuration, the dynamics of the superconducting phase difference across the junction is 
described by the perturbed sine-Gordon equation 



with the boundary conditions 



Vtt - Vxx + sin ip + aipt = 7 + l(t) ! 

7(f) = 7isin(fit) +7 2 sin(2Q£ + 9), (1) 



<p(£) = <p(0) + 2w; <p x (t) = <Px(0) , ( 2 ) 



where x and t are the length along the junction and the time, measured, respectively, in units of the Josephson length 
A j and the inverse plasma frequency (2irf p )~ 1 . The dissipation constant due to the quasiparticle tunneling current, a, 
is of the order of 0.05 in our experiment. Here 7, 71 and 72 are the dc and ac bias current densities, both normalized 
to the junction's critical current density and = w/(27r/ p ) = f/f p is the dimensionless frequency of the ac bias. The 
amplitudes 71 ~ \/P~i and 72 <~ \/p2 correspond to the first and second microwave harmonics, respectively, and 9 is 
the phase shift between them that we controlled in this experiment. 

The effect of a single harmonic microwave (here 72 = 0) on the motion of a fluxon in the junction is presented 
in Fig. 2a. It shows the current-voltage characteristics of the fluxon with no ac drive (open symbols) and with the 
single-harmonic ac drive (solid symbols) having the frequency / = 4.8 GHz. We note that the fluxon characteristics 
in the presence of a monochromatic ac drive of any frequency remains symmetric with respect to the origin, as was 
observed in the earlier experiments 4 ' 5 . 

The effect of the second harmonic in the microwave spectrum is reported in Fig. 2b. The two curves, shown in 
that figure, correspond to two different phase shifts between the drives Pi and Pi- The salient feature of these curves 
is the non-zero voltage at zero dc current through the junction. The sign of this rectified voltage gets reversed by 
changing the phase shift 9 from 7r/2 to 37r/2 (or to —ir/2). The dependence of the rectified voltage at zero bias on the 
phase shift is presented in Fig. 3(a). We observe that the rectified voltage reverses its sign close to 9 = and 9 = ir 
and changes periodically with the phase shift. We also measured the dependence of this effect on power and found 
that, at somewhat lower power levels, the dependence of the rectified voltage on the phase shift may have several 
maxima within period of 27r, see Fig. 3(b). Possibly, at these values of the ac bias the inhomogeneity of the junction 
becomes significant, and the ideal model accounted by Eq. (1) is no longer valid. Numerical simulations of Eq. (1) 
with boundary conditions (2) have been performed with the 4th order Rungc-Kutta method. The average voltage 
drop has been computed and plotted in Figs. 3 (a,b) with the solid line. The following experimental values of the 
systems parameters have been taken: a — 0.05, Q = 0.01, 71 = 72 = 0.1 for the case of Fig. 3(a) and £1 = 0.02, 
71 = 72 = 0.135 for the case of Fig. 3(b). 

The results of the experimental measurements and numerical simulations are in good correspondence with the results 
of the first order (point-particle approximation) soliton perturbation theory 21 . In this approximation the perturbation 
is assumed to be small, so that it does not change the shape of the soliton and only the soliton parameters vary in 
time. Also, this approximation works well in the adiabatic limit, when O <C a. Using this approach, the mean fluxon 
velocity (in the absence of dc bias) can be computed 16 as follows: 

(v) ~7?72Bin(g + fl ), tanflo = r 2 5=-. (3) 

fi[3+(g) 

The voltage drop is proportional to the average fluxon velocity. Thus, the rectified voltage behaves as a sine function 
of the phase shift. Note that #0 - * if — > 0, and this is in accordance with experimental and numerical data, where 
the parameter 9q is small. 

The maximum amplitude of the rectified voltage is rather large - it corresponds to the average normalized fluxon 
velocity of about 0.2. Figure 4 shows the power dependence of the maximum amplitude of the rectified dc voltage. 
For each of the two curves the power level of another harmonic was kept constant (but not referenced). One can see 
that the rectified voltage has a maximum value at some optimum power level. The most important system parameter 
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for the phenomenon is the normal state resistance of the junction. We find that by decreasing the damping in the 
system both the optimal power level and for maximum rectified voltage increase. The dependence of the bi-harmonic 
drive on the phenomenon can be qualitatively undertstood in a simple way. In the limit of the small rf frequency / 
(with respect to the plasma frequency f p ) the kink behaves as a point-like particle and the system reacts to the rf 
drive adiabatically. For any chosen fundamental period of the microwave, the phase shift between two drives makes 
the total current large and positive over short time and then smaller and negative over somewhat longer time. Since 
the current amplitude is large enough to reach the saturation voltage for both current polarities (Swihart velocity 
in the case of fluxon), the mean voltage will be negative in this case. Changing the phase shift simply reverses the 
polarities of the currents and voltage. For higher frequencies, however, the adiabatic approach does not work and 
the phenomenon becomes more complicated. In particular, besides the point particle contribution characterizing the 
adiabatic regime, there is an equally important contribution to the drift current coming from internal oscillations of 
the kink (internal mode mechanism). These oscillations are asymmetric in space and phase locked to the external 
driver and, in presence of damping, they are coupled to the center of mass (translational mode) giving an extra 
contribution to the transport 16 . It would be interesting to separate these two contributions, a task which is not easy 
to perform within our present experimental setting. 

In conclusion, we reported the first experimental observation of the rectified dc voltage, induced by the ratchet-like 
motion of a fluxon in uniform annular Josephson junction driven by bi-harmonic microwave fields. We find that the 
symmetry breaking in the rf field j(t) leads to the unidirectional motion of a fluxon, which is manifested in the non-zero 
rectified voltage across the junction at zero bias current. Our results qualitatively confirm the theoretical predictions 
made in Rcfs. 15, 16. The rectified voltage changes periodically with the phase shift 9. The authors acknowledge 
discussions with P. L. Christiansen, M. Fistul, S. Flach, and M. R. Samuelsen. M. S. acknowledges partial support 
from a MURST-PRIN-2003 Initiative, and from the European grant LOCNET no. HPRN-CT-1999-00163. Y.Z. 
wishes to acknowledge support from the INTAS Young Scientist Fellowship no. 03-55-1799. 
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FIG. 1. Schematic view of an annular Josephson junction with a single fluxon trapped in it. The junction is biased by a dc 
current and an ac bias, induced by applied microwaves. 




' 1 ' 1 ' 

phase shift 

— ■ — ~ ji/2 
• ~ 3rc/2 


' 1 ' 1 1 f 







-60 -40 -20 20 40 60 

dc voltage (nV) 

FIG. 2. (a) Current- voltage characteristics of a single fluxon in annular Josephson junction with no ac drive (open symbols) 
and with single-harmonic ac drive (solid symbols) having the frequency /i = 4.8 GHz. (b) A single-fluxon current-voltage 
characteristics. The fluxon is driven by bi-harmonic ac drive having the frequencies /i = 1.2 GHz and /2 = 2.4 GHz and power 
P2 — +4dB. The phase shift 9 between the two harmonics is fixed to 7r/2 (open symbols) and 37r/2 (solid symbols). 
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FIG. 3. A dependence of the rectified dc voltage measured at zero dc bias current on the phase shift 9 between the two 
harmonics of ac drive, (a) /i = 1.2 GHz and f 2 = 2.4 GHz and power P 2 = +4dB. (b) /i = 2.4 GHz and / 2 = 4.8 GHz, for 
two power levels. Solid lines show results of numerical computation of the voltage for the dimensionless ac bias amplitudes 
E x =Ei = 0.1 (a) and E X =E 2 = 0.135 (b). 
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FIG. 4. A dependence of the maximum amplitude of the rectified dc voltage versus power for the two harmonics of ac drive. 
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